What is claimed is: 



CLAIMS 

1 . A method of correcting data frames for dynamic phase errors in magnetic 
resonance imaging comprising the steps of: 

a) acquiring a plurality of high-resolution data frames, 

b) acquiring with each high-resolution data frame a low-resolution phase reference 
navigator, and 

c) rephasing each high-resolution data frame in A>space with a convolution with 
the Fourier transform of the phase conjugate of the low-resolution phase reference 
navigator. 

2. The method as defined by claim 1 wherein step c) includes the steps of: 

cl) calculating for each data frame a Fourier transform of a phase conjugate of the 
phase reference navigator as a refocusing kernel for the data frame, and 

c2) convolving in A>space the high-resolution data frame with the refocusing 
kernel to phase correct the high-resolution data frame. 

3. The method as defined by claim 2 wherein step c) utilizes a gridding 
reconstruction of each high-resolution data frame. 

4. The method as defined by claim 3 wherein the acquired MRI data, d, is: 

d = GFPm = RPm 

where d is a C R x 1 vector containing the fc-space data, misaA^x 1 vector containing 
the object magnetization in Cartesian coordinates, and P (CN 2 x N 2 ), F (CN 2 x CN 2 ) and 
G (CR x CN 2 ) are matrices respectively representing image-space phase corruption, 
discrete Fourier transform, and resampling from a Cartesian grid onto the &-space 
trajectory. 

5. The method as defined by claim 4 wherein the reconstructed image is a least 
squares estimate: 

m = (P * R * RP) " 1 P * R * d 
= M , P*R*d 
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where * denotes conjugate transpose. 

6. The method as defined by claim 5 wherein the least squares reconstruction is 
estimated as: 

5 m P*R*d. 

7. The method as defined by claim 6 wherein &-space refocusing of each data frame 
includes: 

initialize a zero-filled N x N matrix (which will accumulate the refocused 
spectrum), 

for each interleave: 

(a) Reconstruct the navigator data in /r-space (n x n), 

(b) Zero-pad (a) by a factor of 2 (2/z x 2«), 

(c) Inverse Fourier transform (b), 

(d) Calculate the phase conjugate of (c), 

(e) Fourier transform (d) (to get the refocusing kernel), 

(f) Reconstruct the high-resolution data frame in /:-space at the final image 
resolution (N x TV), 

(g) Convolve (f) with (e) (refocusing step), 

(h) Add (g) to the sum accumulating in step (1), 
inverse Fourier transform the cumulative sum from (h) to get the refocused image. 

8. The method as defined by claim 7 wherein step a) and step b) include acquiring 
information relevant to time-varying data corruption, and ordering of data frames based 

25 on this information thereby promoting smooth, non-periodic modulation of data in the 
spectral domain. 

9. The method as defined by claim 7 and further including before step (c) the step of: 

(bl) weighting each image-space sample in the plurality of data frames by 
30 magnitude of the low-resolution phone reference navigator acquired therewith. 

10. The method as defined by claim 9 wherein 

Qd= QRPm 

where Q is a preconditioning matrix that provides desired image-space weighting. 
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1 1 . The method as defined by claim 10 wherein least-squares solution, m Q , is given 
by: 

m Q =(P*R*Q*QRP)-' P*R*Q*Qd. 

5 

12. The method as defined by claim 1 and further including the step of: 
d) summing the phase-corrected high resolution data frames. 
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